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ABSTRACT 

In the standard paradigm, satellite galaxies are believed to be associated with the population of 
dark matter subhalos. The assumption usually made is that the relationship between satellite galaxies 
and sub-halos is similar to that between central galaxies and host halos. In this paper, we use the 
conditional stellar mass functions of satellite galaxies obtained from a large galaxy group catalogue 
together with models of the subhalo mass functions to explore the consequences of such assumption 
in connection to the stellar mass function of satellite galaxies and the fraction and fate of stripped 
stars from satellites in galaxy groups and clusters of different masses. The majority of the stripped 
stars in massive halos are predicted to end up as intra-cluster stars, and the predicted amounts of 
the intra-cluster component as a function of the velocity dispersion of galaxy system match well the 
observational results obtained by Gonzalez et al. (2007). The fraction of the mass in the stripped 
stars to that remain bound in the central and satellite galaxies is the highest (~ 40% of the total 
stellar mass) in halos with masses Mh ~ lO^"* /i~^Mq. If all these stars end up in the intra-cluster 
component (Max), or a maximum amount of these stars is accreted into the central galaxy (Min), 
then the maximum fraction of the total stars in the whole universe that is in the diffused intra-cluster 
component is ^ 19%, and the minimum is 5%. In the case of 'Max', the stellar mass of the intra- 
cluster component in massive halos with Mh ~ 10^^ /i~"'^Mq is roughly 6 times as large as that of the 
central galaxy. This factor decreases to ~ 2, 1 and 0.1 in halos with Mh ^ 10^**, 10^'^, and 10^^ /i^^Mq, 
respectively. The total amount of stars stripped from satellite galaxies is insufficient to build up the 
central galaxies in halos with masses < lO^^'^ h~^y[Q, and so the quenching of star formation must 
occur in halos with higher masses. In semi-analytical models and simulations that do not resolve 
the diffused component, caution must be exercised when using the observed stellar mass/luminosity 
function of galaxies to constrain the star-formation, feedback and merger processes in dark matter 
halos. 

Subject headings: dark matter - large-scale structure of universe - galaxies: halos 



1. INTRODUCTION 

Recent years have seen a dramatic impetus to link 
galaxies to their dark matter halos. In particular, the 
development of powerful statistical tools such as the halo 
model, the halo occupation distribution (HOD) and the 
conditional luminosity function (CLF), combined with 
the availability of large redshift surveys such as the two- 
Degree Field Galaxy Redshift Survey (2dFGRS; CoUess 
et al. 2001) and the Sloan Digital Sky Survey (SDSS; 
York et al. 2000), have resulted in reliable descriptions 
of how galaxies with different properties are distributed 
over halos of different masses (e.g. Jing et al. 1998; Pea- 
cock & Smith 2000; Berlind & Weinberg 2002; Yang 
et al. 2003; van den Bosch et al. 2003, 2007; Zheng et 
al. 2005; Tinker et al. 2005; Cooray 2006; Brown et 
al. 2008; Cacciato et al. 2008). 

An important aspect of this galaxy-dark matter con- 
nection is that there are two different kinds of galaxies: 
central galaxies, which reside at rest at the center of their 
dark matter halo, and satellite galaxies, which orbit the 
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halo associated with a central galaxy. It is generally be- 
lieved that satellite galaxies themselves reside in dark 
matter subhalos. Prior to being accreted into their cur- 
rent host halo, satellite galaxies were central galaxies, 
and their associated subhalos were host halos (through- 
out this paper we use the term 'host halo' to refer to a 
dark matter halo which is not a subhalo). This implies 
a direct link between the occupation statistics of dark 
matter halos, and their merger/ accretion histories. 

While orbiting the host halo, subhalos and their associ- 
ated satellite galaxies are subjected to dynamical friction 
which causes the substructure to loose its momentum 
and to sink towards the center of the host halo. In the 
mean time, the substructure is subjected to tidal forces 
which try to dissolve it. In particular, tidal heating and 
stripping cause the system to loose mass, and may even 
result in the complete disruption of a subhalo and its 
satellite galaxy. There are thus three possible fates for 
the stars in satellite galaxies: (i) they remain bound to 
a surviving satellite galaxy, (ii) they are accreted by the 
central galaxy, or (iii) they are stripped from the satellite 
galaxy, giving rise to a stellar halo, which are stars that 
orbit the host halo but that are not gravitationally bound 
to any particular galaxy. Throughout this paper we will 
refer to the stars that belong to (ii) and (iii) combined as 
the 'non-surviving population'. Obviously, the satellites 
that are disrupted contribute their entire stellar content 
to the central galaxy or stellar halo, while the survived 
satellite galaxies may also contribute significantly due to 
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tidal stripping. 

The idea that sateUite galaxies are stripped and/or dis- 
rupted is a standard prediction of hierarchical models of 
structure formation. Numerous studies have addressed 
how this phenomenon gives rise to stellar halos in sys- 
tems ranging from spiral galaxies like our own Milky 
Way (e.g. Searle & Zinn 1978; Jonston et al. 1996; 
2001; Robertson et al. 2005; Font et al. 2006) to rich 
galaxy clusters (e.g. Gallagher & Ostriker 1972; Mer- 
ritt 1983; Mihos 2004; Willman et al. 2004; Lin & Mohr 
2004; Conroy et al. 2007; PurceU et al. 2007; 2008; Hen- 
riques et al. 2008). There is also ample observational 
support for the existence of stellar halos formed out of 
disrupted satellite galaxies. In particular, in recent years 
it has become clear that the stellar halo of the Milky Way 
reveals a large amount of substructure in the form of stel- 
lar streams (Helmi et al. 1999; Yanny et al. 2003; Bell 
et al. 2008), In some cases these streams can be unam- 
biguously associated with their original stellar structure 
(Ibata, Gilmore & Irwin 1994; Odenkirchen et al. 2002). 
Similar streams have also been detected in our neigh- 
bor galaxy M31 (e.g. Ferguson et al. 2002). Unfortu- 
nately, due to their extremely low surface brightnesses, 
it is very difficult to detect stellar halos in more distant 
galaxies (but see Zibetti et al. 2004). Consequently, our 
knowledge of stellar halos around individual galaxies is 
fairly limited. However, groups and clusters of galaxies 
also contain a significant stellar halo component, which is 
usually referred to as 'intra-cluster stars' (ICS)^. In fact, 
data indicates that the fraction of stars associated with 
such an ICS component increases with increasing halo 
mass (Gonzalez, Zabludoff & Zaritsky 2007): in massive 
clusters the mass of the stellar halo can be as large as ten 
times the stellar mass of the central (brightest) cluster 
galaxy (Gonzalez, Zabludoff & Zaritsky 2005; Seigar et 
al. 2007). 

With the use of high-resolution numerical simulations, 
it has recently become possible to accurately determine 
the properties (mass function, spatial distribution, or- 
bits, density profiles, spins) of the population of dark 
matter subhalos (e.g. Gao et al. 2004; De Lucia et 
al. 2004; Tormen et al. 2004; van den Bosch et al. 2005; 
Weller et al. 2005; Diemand et al. 2007; Giocoh et 
al. 2008). If subhalos are indeed associated with satellite 
galaxies, these properties should be related to the lumi- 
nosity function, spatial distribution, orbits, and struc- 
tural properties of satellite galaxies. Unfortunately, the 
exact link between satellite galaxies and dark matter 
subhalos is not trivial. Since the stellar component of 
a satellite is expected to be more tightly bound than 
its surrounding dark matter (due to dissipation during 
the formation process), the dark matter is more easily 
stripped, thus causing the ratio between stellar mass and 
total mass to increase with time. Consequently, it is to 
be expected that the occupation statistics of subhalos as 
a function of their (current) mass are different from those 
of host halos. However, since subhalos were host halos 
before being accreted, it seems likely that the occupa- 
tion statistics of subhalos as a function of their mass at 
the time of accretion are identical to those of host ha- 
los at that time. Indeed, it has been shown that models 

^ Throughout this paper we will use the terms 'ICS' and 'stellar 
halo' without distinction. 




Fig. 1. — The solid line shows the 'uncvolvcd' subhalo mass func- 
tion fraction predicted by Giocoli et al. (2008) from N-body simu- 
lations. The dotted and dashed lines are the predicted mass func- 
tions of sub^-subhalo and sub^ -subhalos, respectively (see Eq.[[2]). 

based on this ansatz are extremely successful in explain- 
ing the correlation functions and luminosity functions 
of galaxies at different redshifts (e.g. Vale & Ostriker 
2004, 2006; Conroy, Wechsler & Kravtsov 2006). How- 
ever, since satellite galaxies only make up a small frac- 
tion ( < 30%) of the total galaxy population (e.g. van 
den Bosch et al. 2007, 2008; Tinker et al. 2007; Yang et 
al. 2008a), this consistency cannot be considered a sen- 
sitive test of the subhalo - satellite connection. This is 
also apparent from the fact that models that link satel- 
lite properties to the current subhalo mass can also fit the 
data remarkably well (see e.g., Mandelbaum et al. 2006; 
Kim et al. 2008). 

In this paper, we use the conditional stellar mass func- 
tion of satellite galaxies and the relation between stellar 
mass and halo mass for central galaxies, together with 
the subhalo mass functions obtained from recent numer- 
ical simulations, to determine the relation between satel- 
lite galaxies and dark matter subhalos. In particular, we 
investigate what fraction of satellite galaxies survives, 
what fraction is accreted into the central galaxy, and 
what fraction is tidally disrupted. Our approach is em- 
pirical, because the stellar mass functions and the cen- 
tral stellar mass - halo mass relation are obtained from 
a large SDSS galaxy group catalogue. The structure 
of this paper is organized as follows. Section [2] gives 
a brief description of the data and the model used in 
this paper. In Section [3] we present our predictions of 
the conditional stellar mass functions for satellite galax- 
ies. In Section [1] we discuss the possible fates of the 
stars that are stripped from the satellite galaxies. Fi- 
nally, we summarize our results in Section (5] Through- 
out this paper we adopt a ACDM cosmology with pa- 
rameters that are consistent with the three- year data re- 
lease of the WMAP mission (hereafter WMAP3 cosmol- 
ogy): = 0.238, Qa = 0.762, f^b = 0.042, n = 0.951, 
h = iJo/(100 kms~iMpc"^) = 0.73 and as = 0.75 
(Spergel et al. 2007). 

2. DATA AND ANALYSIS 
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Fig. 2. — The conditional stellar mass functions for satellite galaxies in halos of different masses as indicated in each panel. The circles 
are the measurements from the SDSS DR4 group catalogue. The error-bars represent scatter among four different samples as described in 
the text. Dotted, dashed and long-dashed lines in each panel show the predicted CSMFs using the unevolved subhalo, sub^-subhalo and 
sub-^-subhalo mass functions, respectively. Here we assume that the relation between stellar mass and halo mass is the same for satellite 
galaxies and central galaxies. The solid lines are the combinations of the three components. Note that we did not take into account the 
contribution to the satellite population, which is negligible, from higher order subhalos, e.g., sub^-subhalos and sub'^-subhalos etc. 



The data used in this paper is based on the SDSS 
galaxy group catalogues of Yang et al. (2007; hereafter 
Y07). These catalogues are constructed by applying the 
halo-based group finder developed by Yang et al. (2005) 
to the New York University Value-Added Galaxy Cat- 
alogue (NYU-VAGC; see Blanton et al. 2005), which is 
based on SDSS Data Release 4 (Adelman-McCarthy et 
al. 2006). Detailed tests have shown that this group 
finder is very successful in associating galaxies accord- 
ing to their common dark matter halos, and that the 
halo masses that have been assigned to the groups are 
reliable. 

In Yang et al. (2008a,b) we have used these group 
catalogues to determine the conditional luminosity func- 
tions (CLF) and the conditional stellar mass functions 
(CSMF) separately for central and satellite galaxies. 
Central galaxies are defined as the most massive group 
members in terms of their stellar mass, and satellites are 
the group members that are not centrals. The CSMF of 
satellite galaxies, 'i>s(-^*|-Mh), gives the average number 
of satellites with stellar masses in the range M» ± dM* 
that reside in a host halo of mass Mh. The open cir- 
cles in Fig. [2] show the $s(-M*|-Mh) obtained by Yang et 
al. (2008b, hereafter YMB08) for halos in four different 
mass bins. These CSMFs are the averages obtained from 
two group samples (Samples II and III) and from using 
two different halo mass estimators (see YMB08 for de- 
tails). The errorbars in Fig. [5] reffect the scatter among 
these four results, which in general are much larger than 



the statistical errors obtained using bootstrap samples. 

In order to link the satellites to the subhalo popula- 
tion, we also need the subhalo mass function (SHMF). 
In fact, since it is expected that the properties of satellite 
galaxies are linked to the mass of their subhalos at their 
time of accretion (see discussion in SJT]), we need the so- 
called "unevolved" SHMF, n^n{'ms\M\-^) , which gives the 
number of subhalos that have been accreted by the main 
progenitor of a halo of present day mass Mh, as a func- 
tion of their mass nis at the time of accretion. This un- 
evolved SHMF should not be confused with the "evolved" 
SHMF, ncv{m'^\Mii), which gives the number of subha- 
los with present-day masses in the range ± dm^ that 
reside, at present, in a host halo of mass Mh. The un- 
evolved SHMF evolves into the evolved SHMF due to 
the combined effect of dynamical friction, tidal stripping 
and tidal heating, which causes some subhalos to dis- 
solve, and others to loose mass (see van den Bosch et 
al. 2005 for details). 

Using high-resolution numerical simulations, Giocoli et 
al. (2008) found that the unevolved SHMF is accurately 
described by 



?lun,o(n^s|Mi 



0.176 



exp 



-12.27 

VMh 



(1) 

Note, however, that this SHMF only includes the pop- 
ulation of subhalos that are directly accreted onto the 
main progenitor of the host halo (hence the subscript 
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Fig. 3. — Same as Fig. [2] but here we show the predicted total CSMFs using the unevolved (solid lines) SHMF and the best fit results 
(long-dashed lines) of the model described in the text, respectively. 



'0'). This neglects the possibility that subhalos may 
themselves contain 'sub-subhalos', and sub-subhalos may 
contain 'sub-sub-subhalos' and so on. For convenience, 
we will use the notation 'sub'-subhalo' {i — 1,2,3, 4...) to 
refer to the i^^ level of subhalos. In modeling the satel- 
lite population associated with the subhalos, neglecting 
the 'sub'-subhalo' populations may result in an underes- 
timate of the number of satellite galaxies. This is espe- 
cially true for massive systems, where the merging pro- 
genitors may already contain relatively massive satellite 
galaxies. 

Fortunately, because of the self-similar behavior of the 
unevolved SHMF (Eq. [T]) , we can calculate the unevolved 
SHMF including the sub'-subhalo populations. Assum- 
ing that Eq. ^ also apphes to subhalos, the SHMF of 
sub'-subhalos can be written as 



nun,iims,t\Mh) 



i_i(ms|Mh)dms 



(2) 

Fig. [1] shows the predictions for the unevolved subhalo, 
sub^-subhalo and sub^-subhalo mass functions, thus ob- 
tained. As one can see, the contributions from sub-- 
subhalos and sub^-subhalos to the total SHMF can be 
larger than that from the subhalos at log(TOs/Mh) < —3 
and < —6, respectively. They can thus contribute a sig- 
nificant fraction of (small) satellite galaxies. 

The final ingredient for our modeling is the conditional 
probability distribution, VdM^lMy^), that a halo of mass 
Afh hosts a central galaxy with stellar mass Using 
our SDSS galaxy group catalogues, YMB08 found that 
7'c(Af* I A/h) is well described by a log-normal distribution 



whose median is given by 

{M,)iAh) = Mo 



'(l + Mh/Mi)/3- 
Here Mi is a characteristic halo mass so that cx 
M^'^^ for A/h < Afi and Af, cx M^ for Afh > A/i. 
The parameters obtained from the SDSS groups are: 
logAfo = 10.306, logAfi = 11.040, a = 0.315, and (3 = 
4.543, where M^ is in units of R-'^Mq and Mi in /i'^Mq 
(see YMB08 for details). The width of Vc{M^\My,) is 
found to be roughly independent of halo mass (see also 
More et al. 2008), with a dispersion cr(logAf*) — 0.173. 

3. THE DISRUPTION OF SATELLITE GALAXIES 

The prediction of the CSMF of satellite galaxies can 
be written as 



*,(Af,|Afh) 



7's(Af*|ms)nun('ns|Afh)dms , (4) 



where Vs{M^,\ms) is the probability that a subhalo of 
mass rus at the time of accretion hosts a present-day 
satellite galaxy with stellar mass Af* , and 



nun(n^s|Afh) 



?^un,^(n^■s|A^h) 



(5) 



is the unevolved SHMF, including all sub*-subhalos up to 
level iVinax- We start by assuming that the stellar mass - 
halo mass relation does not evolve with redshift, so that 
'Ps[M^\ms) — Vc{M^.\M\i — nis). In this case, we denote 
the model prediction of <I>s by $un, and so 



$u„(Af*|A//h) 



Vc (A/* I TOs ) Uunims I A'/h ) dms . (6) 
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Fig. 4. — The best-fit (survival) fraction of subhalos that can 
host sateUite galaxies in the same way as halos can host central 
galaxies at present time 2 = 0. Here we assume that the subhalos 
in host halos of different masses have the same survival fraction as 
a function of ms/Afjj. As shown in Fig. [2]the best fit predictions 
of this model are remarkably good. 

The solid lines in each of the panels of Fig. [2] show 
the CSMFs thus obtained with Eq. ^ and using 
-^max — 2. The dotted, dashed and long-dashed lines 
show the contributions due to subhalos, sub^-subhalos 
and sub^-subhalos, respectively, obtained by replacing 
"■un(ws|Mh) in Eq. ^ with nun,i(™s|^h) for i = 0, 1 
and 2. Note that for massive halos (A4 > 10^"* /i^^Mq) 
the contribution of sub^-subhalos is comparable to that 
of the subhalos, indicating that their contribution can- 
not be ignored. In low mass halos (Mh < 10^^ h~^MQ), 
on the other hand, ignoring the sub*-subhalos for i > 1 
makes little difference. The satellite population con- 
tributed by the sub^-subhalos is roughly 10% in massive 
halos with Mh ~ lO" H'^Mq but drops to < 1% in 
halos with Mh ~ 10^^ /i"^Mq. We found that the over- 
all contribution by sub'-subhalos with i > 2 is always 
less than a few percent and can be safely neglected. In 
what follows, we therefore adopt iVmax = 2 throughout, 
and unless stated otherwise, we refer to subhalos, sub^- 
subhalos and sub^-subhalos collectively as 'subhalos'. 

Comparing the model predictions with the observa- 
tional data (open-circles with error-bars) shown in Fig. 
[21 we see that the predicted CSMFs for satellite galax- 
ies are roughly consistent with the data for massive 
groups. However, the model significantly over-predicts 
the CSMFs for low-mass halos. There are two different 
possible explanations for this discrepancy. First, halos 
of a given mass at high redshift may contain different 
amounts of stars from halos at z = with the same mass, 
so that the stellar mass - halo mass relation for subhalos 
is different from that of present-day central galaxies. In 
order to explain the discrepancy, the stellar mass frac- 
tions of halos of a given mass then have to be lower at 
higher redshifts. Furthermore, since the extent of the 
discrepancy depends on halo mass, the amplitude of the 
redshift dependence has to be different for halos of dif- 
ferent masses. This solution is not very likely, since it 
is rather contrived to assume that the star formation ef- 
ficiency in progenitor halos depends on the mass of the 



halo in which it will end up in the future. Nevertheless, 
we acknowledge that the average relation between halo 
mass and stellar mass may well evolve with redshift (see 
e.g., Conroy & Wechsler 2008 for empirical constraints 
in support of such evolution). However, given the typ- 
ical accretion times for subhalos, we believe that this 
will not have a strong impact on our results. This is 
also supported by the work of Purcell et al. (2007), who 
have shown that the predictions for the stellar halo mass 
fractions are extremely robust to changes in the star for- 
mation histories of the galaxies. The second, more likely, 
possibility is that the subhalos and their satellite galax- 
ies experience mass loss and/or disruption due to the 
combined effect of dynamical friction and tidal forces. 

In this paper we focus on the second possibility. In 
this case, a satellite galaxy can either be completely dis- 
rupted, hence does not contribute to the satellite popu- 
lation, or experience mass loss but survives as a satellite 
of lower mass. Both of these effects can change the pre- 
dicted stellar mass function of satellite galaxies relative 
to <I>un- As an illustration, let us consider a simple model 
in which a satellite is either completely disrupted or re- 
mains intact. This assumption is consistent with the 
simulation results that a satellite is quickly disrupted af- 
ter it has lost a significant amount of mass (Moore et al. 
1999). It is also valid if the disruption of a satellite is 
due to merger into the central galaxy. As shown in van 
den Bosch (2005) and Giocoli et al. (2008), on average 
the instantaneous mass loss rate of a subhalo depends on 
the ratio between the instantaneous subhalo mass and 
the host halo mass at the time in question. Here we ig- 
nore such details. Instead we consider a simple model 
where the survivor fraction, /, is a function of the ra- 
tio between the subhalo mass at accretion, rris, and the 
host halo mass at the present time, Mh. Motivated by 
the roughly self-similar behavior of the subhalo popu- 
lation, we assume that /(mg/Mh) is universal in halos 
of different masses. For a given /(mg/Mh) the CSMF 
is given by Eq. (|4]), but with ?iun(TOs|-^h) replaced by 
/(mg/Mh) nun(ws|Mh). We model f{x) using a polyno- 
mial form, f{x) = a + bx + cx^ + dx^ , and determine the 
free parameters (a, b, c, d) by fitting the model to nine 
CSMFs that cover the mass range 12.0 < logMh < 14.7, 
each with a 0.3 dex bin width in halo mass. The CSMFs 
corresponding to the best-fit model are shown as the 
long-dashed lines in Fig. [3l Note that this simple model 
fits the data surprisingly well, supporting the assumption 
that the survivor fraction depends only on mg and Mh 
via their ratio. The best-fit /(ms/Mh) is shown in Fig.lH 
which shows that the survivor fraction decreases mono- 
tonically with subhalo mass, from ^ 1 for subhalos with 
nis ~ 10~'^-^Mh to ~ 0.1 for subhalos with nig ~ 0.6Afh 
(recall that mg refers to the subhalo mass at the time of 
accretion). Thus, if the stellar mass - halo mass relation 
does not evolve with redshift, then more massive (rela- 
tive to their host) subhalos and their associated satellite 
galaxies are predicted to have a smaller survival probabil- 
ity. This is consistent with a picture in which dynamical 
friction is responsible for transporting satellite galaxies 
to the inner regions of their host halos, where they are 
more likely to be disrupted by tidal forces or to merge 
with the central galaxy. 
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Fig. 5. — Left panel: the long-dashed line shows the predicted average total stellar mass that the satellite galaxies (in unevolved subhalos) 
can bring into the host hales. The dashed line shows the average total stellar mass remain in the satellite galaxies which can be observed 
today. The dot-dashed line corresponds to the average total stellar mass that has been stripped or disrupted from the satellite galaxies. 
For comparison, the solid line shows the predicted average stellar mass of the central galaxy. Right panel: the ratios of the stripped or 
disrupted stellar mass from satellite galaxies, Af,_ns, to that of the central galaxy, Mt,c, (solid line) and to the total stellar mass in central 
and satellite galaxies, M»^c + Af»^s, (dashed line). 



4. MERGING, TIDAL DISRUPTION AND 
QUENCHING 

We can use the stellar mass - halo mass relation to 
predict the following stellar mass components for halos 
of a given mass: the total stellar mass that the satellite 
galaxies bring into the host halo, which is given by 

/•CO 

M,,u„(Afh)= / dM,Af,$u„(Af*|Mh), (7) 
Jo 

and the total stellar mass in surviving satellite galaxies, 

/•OO 

M,.,(Mh)= / dA/,M,$,(M,|Mh). (8) 
Jo 

Note that both M*,u„(Mh) and M*^s(Mh) can be ob- 
tained without assuming whether a satellite is completely 
disrupted or only partly stripped. The difference between 
these two stellar masses, Af*^ns(Afh) = Af*,un(A'/h) — 
A<f*,s(A'/h), is the total stellar mass of the non-survivors, 
which consists of the stars in satellite galaxies that are 
completely disrupted and those that are stripped from 
the satellite galaxies. In other words, M*_ns is the sum of 
the stellar mass of the stellar halo (i.e., the ICS) plus the 
stellar mass of the central galaxy that has been accreted 
from satellite galaxies. 

The left-hand panel of Fig. [5] shows the predictions for 
A^»,un (long-dashed line), M*,s (short-dashed hue), and 
A<r*_ns (dot-dashed line), all as functions of halo mass. 
For comparison, the solid line shows the stellar mass of 
central galaxies as a function of halo mass [i.e. the stellar 
mass -halo mass relation given by Eq. Note that, in 
massive halos with Mh > 10^'^ h~'^yiQ, the stellar mass 
in the non-surviving component is larger than the stellar 
mass of the central galaxy. 

The dashed line in the right-hand panel of Fig.[5]shows 
A<r*,ns/(A^*,c + Af*,s), the ratio between the stellar mass 
of the non-surviving component and the combined stel- 
lar mass of the central galaxy and the surviving satel- 
lites. In terms of Af*_ns/ (Af*.c + Af^.s), halos with Afh ^ 
H-^Mq have the highest ratio, - 0.6, i.e. ^ 40% 



of the total stellar mass is in the non-surviving compo- 
nent. In less massive halos, the ratio drops rapidly due 
to the fact that the stellar mass - halo mass relation is 
extremely steep at the low mass end, so that low mass 
subhalos contain very few stars even at their time of ac- 
cretion. At the massive end the ratio also decreases due 
to two effects. First of ah, for My, > 10^^ /i-^Mg the 
stellar mass fractions decline with increasing halo mass 
(i.e., the slope of the stellar mass - halo mass relation 
is less than unity), so that more massive subhalos con- 
tribute fewer stars per unit dark matter mass. Secondly, 
many of the fainter satellite galaxies, which reside in less 
massive subhalos, can survive since dynamical friction is 
not effective for subhalos with small nis/Mh- 

The ratio of stars that are in the non-surviving com- 
ponent in the entire Universe to those that are locked 
up in either central or satellite galaxies can be estimated 
using: 

^ /o'"Af.,ns(Afh)n(A//h)dA/h 

" /o°° [M.AMh) + Af,,,(Afh)] n{A'h)dAh ' ^ ' 

where ri(Afh) is the halo mass function (e.g., Sheth, Mo 
& Tormen 2001; Warren et al. 2006). We find that 
Rns ~ 23%. If we assume that all of the stars in the non- 
surviving component are turned into ICS, then ~ 19%, 
which is the upper-limit, of the stellar mass in the uni- 
verse is in ICS. On the other hand, if we assume that 
maximum of these stars are accreted by the central galax- 
ies, then ~ 5%, which is the lower-limit, of the stellar 
mass in the universe is in ICS. 

4.1. The fate of stripped and dispersed stars 

We now investigate the fate of the stars that are 
stripped from satellite galaxies. In particular, we exam- 
ine what fraction is actually disrupted by tidal forces, 
thus giving rise to a stellar halo (i.e. the ICS), and 
what fraction is accreted by the central galaxy. For 
this purpose, we show as the solid line in the right-hand 
panel of Fig.[5l the ratio between the stellar mass of the 
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Fig. 6. — The fraction of the total stellar masses, which is 
contained in the central galaxy and intra-cluster stars as a function 
of cluster velocity dispersion. Data with errorbars are obtained 
from Gonzalez et al. (2007), measured within rsoo for the 23 groups 
and clusters in their sample. The solid and dashed lines are our 
model predictions for the two extreme cases. See text for details. 

non-surviving component, Af^ „s, and that of the central 
galaxy M» c- As one can see, in small mass halos with 
Mh < 10^^h~^MQ, the mass of the stripped stars is much 
smaller than that of the central galaxy, indicating that 
the latter cannot have grown substantially due to the ac- 
cretion of satellites. Rather, central galaxies in low mass 
halos must have grown predominantly via star formation. 
It also means that if all the mass of the stripped stars 
ends up as a stellar halo, the mass of that stellar halo 
can only be a small fraction of the mass of the central 
galaxy. In halos with Mh ^ 10^^ H'^Mq, on the other 
hand, M*^„s S> M^^c, and a very significant fraction of 
the stellar mass of the central galaxy may consist of ac- 
creted stars stripped from satellites (but does not have 
to). It is also clear that we can exclude the possibil- 
ity that all stripped stars are accreted into the central 
galaxy, as this would imply stellar masses for the bright- 
est cluster galaxies that are about 8 times higher than 
observed. Rather, the stripped stars must have given rise 
to a substantial stellar halo. 

The solid dots with errorbars in Fig. [S] show the ob- 
served fractions of the total stellar mass present in groups 
and clusters that is contained in the central galaxy and 
the ICS as a function of the line-of-sight velocity dis- 
persion of the member galaxies. These data are taken 
from Gonzalez, Zaritsky & Zabludoff (2007), under the 
assumption that satellite galaxies, central galaxies and 
the ICS all have the same mass-to-light ratios (in the 
i-band). We can compare the data with two extreme 
predictions of our model, by converting halo mass into 
a line-of-sight velocity dispersion using equation (6) in 
Yang et al. (2007). The solid lines, labeled 'Max', corre- 
sponds to the prediction under the assumption that all 
stars of the non-surviving component are added to the 
stellar halo (i.e., no stars stripped from satellite galaxies 
are accreted by the central galaxy). Clearly, this corre- 
sponds to the maximum amount of ICS possible. The 
dashed lines, labeled 'Min', corresponds to the minimum 
amount of ICS, i.e. the stellar mass in the ICS is as- 



sumed to be max(M*^ns — Af*_c,0). Note that the data 
of Gonzalez et al. (2007) are obtained within rsoo, the 
radius within which the cluster mass density exceeds the 
critical value by a factor of 500. In our model predic- 
tion, the total mass of satellite galaxies within rsoo is 
estimated by assuming that the distribution of satellite 
galaxies follow the NFW (Navarro, Frenk & White 1997) 
profile with concentration appropriate for the halo mass 
in question. For the distribution of ICS, we consider two 
cases. Case I assumes that the ICS have the same dis- 
tribution as the satellite galaxies, and the corresponding 
results are shown in Fig. [6] as the two thin lines. Case 
II assumes that all ICS are distributed within rsoo, and 
the corresponding results are shown in Fig. [S] by the two 
thick lines. It is reassuring that, for a given assumption 
of the ICS distribution, the two extreme models, 'Max' 
and 'Min', give quite similar results. The prediction of 
Case I is in good agreement with the data over the en- 
tire range of masses probed, while Case II over-predicts 
the ratio moderately. It is unclear whether this moder- 
ate discrepancy necessarily implies that the distribution 
of the ICS extends beyond r5oo, because there are other 
factors that can cause such discrepancy. For example, 
stars in the central galaxy and in the ICS component 
may, on average, have a higher stellar mass-to-light ratio 
than those in the surviving satellites, so that the ratio 
(M*^c + -Af*jcs)/-^*, total is underestimated in the data 
due to the assumption of the same stellar mass-to-light 
ratio for all components. In addition, there are also un- 
certainties in the estimates of the velocity dispersions, 
(7, particularly for poor systems. Given these uncertain- 
ties, we consider the overall agreement between the data 
and the model prediction remarkable. Clearly, better ob- 
servational data are required in order to distinguish the 
different models considered here. 

The idea that a significant fraction of the stars stripped 
from satellite galaxies end up as ICS is not only consis- 
tent (in fact, required) by the data, but has also been 
found in hydrodynamical simulations of galaxy clusters 
(e.g. Napohtano et al. 2003; Murante et al. 2004, 2007; 
Willman et al. 2004; Sommer-Larsen, Romeo & Portinari 
2005; Rudnick, Mihos & McBride 2006). In addition, as 
shown by Monaco et al. (2006) and Conroy et al. (2007), 
the creation of a significant ICS component is also re- 
quired in order to reconcile the low rate at which mas- 
sive galaxies have grown since z = 1, with the relatively 
high rate at which their host halos (i.e., the clusters) 
have grown in mass (see also Brown et al. 2008). Finally, 
Kang & van den Bosch (2008) argue in favor of the cre- 
ation of stellar halos in order to prevent central galaxies 
from accreting too many blue and/or gas rich satellites, 
which would cause central galaxies to be too blue. 

Recently Purcell et al. (2007) used the analytic model 
for subhalo infall and evolution of Zentner et al. (2005), 
combined with empirical constraints on the stellar mass 
fractions of accreted subhalos, to predict the diffuse stel- 
lar mass fractions of dark matter halos. They predict 
that the average stellar mass fraction in diffuse, intra- 
halo light increases strongly from ~ 0.005 for small 
galaxy halos (~ 10^^ h~^MQ) to ~ 0.2 for poor groups 
(~ lO^'^ft.~"'^M0), after which the trend with mass flattens 
considerably. A comparison with the right-hand panel of 
Fig. [5] shows that this is in remarkable agreement with 
our predictions, which have been deduced using a com- 
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Fig. 7. — The ratios of the stripped or disrupted stellar mass 
from satellite galaxies together with the stellar mass of central 
galaxies at the quench point, M,,ns + Aft^Oi to that of the central 
galaxy, Mt,c- See text for details. 

pletely different approach. In addition to this, Purcell et 
ah (2007) also compared their model predictions to the 
observational data of Gonzales et al. (2007), similar to 
our comparison in Fig. [51 and reached a similar conclu- 
sion as ours, namely that the vast majority of subhalo 
stars must be deposited into a diffuse halo component. 
Similar results have also been obtained by Henriques et 
al. (2008) using a semi-analytical model for galaxy forma- 
tion that includes a simple treatment for the disruption 
of satellite galaxies. It is reassuring that different ap- 
proaches yield results that are in such good agreement. 

4.2. A Lower Limit on the Halo Mass for Quenching 
Star Formation 

In hierarchical models of structure formation, massive 
halos are built up by the mergers of smaller ones. Over 
the years, it has become clear that a successful model 
for galaxy formation in such a hierarchical framework 
requires a mechanism that can quench star formation in 
massive galaxies. Currently, the most favored quenching 
mechanism is feedback from an active galactic nucleus 
( AGN) , which is often assumed to operate above a given 
halo mass, Afq (e.g. Croton et al. 2006; Cattaneo et 
al. 2006, and references therein). We can use the results 
presented above to put a lower limit on this quenching 
mass Mq, as follows. 

In general, the stars in a central galaxy may come 
from two different channels: (i) in situ formation, and 
(ii) accretion from satellite galaxies (those that are part 
of the non-surviving component). Consider a model in 
which star formation is abruptly quenched in halos with 
Mh > Mq. In such quenched halos, the total amount 
of stars from channel (i) is at most equal to JVh^, the 
stellar mass contained in the central galaxy of a halo 
with mass Mq. It may be smaller than M,^Oi because 
part of this mass may be due to accretion rather than in 
situ formation. The total amount of stars from channel 
(ii) is at most M^^ns, because part of the disrupted mass 
may end up in the ICS. In Fig. [7l we show the ratio 
(M*_ns + M,^o)/M,_c, as a function of halo mass. Differ- 
ent lines show the results for different values of Afq [— 



(solid line) , lO" ^, 10i2 "^ lO^^-^ 1012-4/^-11^0]. As one 
can see, if Mq < IO^'^-'^H^^Mq, the total stellar mass from 
the two channels, M<,^ns -|- M,^Oj is insufficient to account 
for the stellar masses of the central galaxies in halos with 
Mh - 1012-5 H-^Mq. However, for Mq > 10^2-4 H-^Mq 
the ratio (Af^.ns 4- Af*^o)/M,_c > 1 for all halo masses. 
This implies that if a quenching mass Afq exists, it must 
be larger than or equal to lO^^ "^ H'^Mq. 

5. DISCUSSION AND SUMMARY 

Using the conditional stellar mass functions for satel- 
lite galaxies obtained by Yang et al. (2008b) and the 
sub-halo mass functions given by Giocoli et al. (2008), 
we study the connection between subhalos and satellite 
galaxies. Assuming that at the time of accretion satel- 
lite galaxies are associated with subhalos according to 
the same stellar mass - halo mass relation as present-day 
central galaxies with halos, we predict the stellar mass 
function of satellite galaxies and compare our model pre- 
dictions with observations. Our main results can be sum- 
marized as follows; 

• Assuming that the stellar masses of satellite galax- 
ies do not evolve after they are accreted into the 
host halos, we find that the model over-predicts 
the population of satellite galaxies, especially in 
low mass halos. One solution, albeit unlikely, is 
that the stellar mass fractions of haloes at higher 
redshifts are lower than at the present. The other, 
which is considered in the present paper, is that a 
significant fraction of satellite galaxies have been 
stripped of their stars or even totally disrupted. 

• Assuming that the amount of disruption of satel- 
lite galaxies is a function of the ratio between the 
mass of the subhalo and that of the host halo, we 
find that the surviving fraction can be described by 
f{x = log(ms/Af,0) = -0.041 -0.820a; -0.422a;2- 
0.078a;'^. The decrease of / with increasing nis/Mh 
is consistent with the idea that dynamical friction 
brings subhalos and satellite galaxies to the inner 
part of the host halo where stripping and disrup- 
tion are more efficient. 

• The majority of the stars stripped from satellites 
in massive halos are predicted to end up as intra- 
cluster stars, and the predicted amounts of the 
intra-cluster component as a function of the ve- 
locity dispersion of galaxy system match well the 
observational results obtained by Gonzalez et al. 
(2007). The fraction of the mass in the stripped 
stars to that in the surviving central and satellite 
galaxies is predicted to be the highest (~ 40% of 
the total) in halos with masses ~ lO^^ H^^Mq. If 
all these stripped stars end up in the intra-cluster 
component (Max), or maximum of them are ac- 
creted into the central galaxy (Min), then we can 
predict that a maximum ~ 19% and a minimum 
^ 5% of the total stars in the whole universe are 
in terms of the diffused intra-cluster component. 

• Stars stripped from satellite galaxies are not suf- 
ficient to build up the central galaxies in halos 
with masses < 10^2-5 H^^Mq, and so star forma- 
tion should not be quenched in halos with masses 
up to at least ~ 10i2-4 H-'^Mq. 
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It should be pointed out once again that our results 
are based on the assumption that subhalos have the same 
stellar mass - halo mass relation as the host halos at red- 
shift z — 0. If the total amount of stars that can form in 
a halo at higher redshift is larger than that in a halo of 
the same mass at lower redshift (e.g., Cooray 2005), the 
predicted disrupted fraction of satellite galaxies would 
be larger, and so would be the predicted amount of ICS. 
On the other hand, if the total amount of stars that can 
form in high-z halos is lower than that in its z = coun- 
terparts, the predicted ICS fraction would be lower. It 
is interesting that the observed ICS fraction is matched 
well with the assumption that the total amount of stars 
that can form in a halo of a given mass is independent 
of redshift. This assumption is also consistent with the 
results obtained by Wang et al. (2006) based on a semi- 
analytical model of structure formation. However, cur- 
rent semi-analytical models of galaxy formation ignore 
the existence of the ICS component, so that the stars 
that formed in subhalos are assumed either to remain in 
satellite galaxies or to merge into central galaxies. Con- 
sequently, such models either predict too high a mass for 
the central galaxies in massive halos or over-predict the 
number of satellite galaxies in group-sized halos. 

As described above, the mass fraction of stars in the 
diffuse ICS component is much larger than that in the 
central galaxies for halos with masses > 10^^h~^MQ, and 
even in the whole universe, this stellar mass component 



Adelman-McCarthy J.K., et al., 2006, ApJS, 162, 38 

Bell E.F., et al., 2008, ApJ, 680, 295 

Berlind A.A., Weinberg D.H., 2002, ApJ, 575, 587 

Blanton, M.R. et al., 2005, AJ, 129, 2562 

Brown M.J.I., et al., 2008, preprint ( arXiv:0804.2293 ) 

Cacciato M., van den Bosch F.C., More S. Li R., Mo H.J., Yang 

X., 2008, MNRAS, submitted (arXiv:0807.4932 ) 
Cattaneo A., Dekel A., Devriendt J., Guiderdoni B., Blaizot J., 

2006, MNRAS, 370, 1651 
Colless M., et al., 2001, MNRAS, 328, 1039 

Conroy C, Wechsler R.H., Kravtsov A.V., 2006, ApJ, 647, 201 

Conroy C, Wechsler R.H., Kravtsov A.V., 2007, ApJ, 668, 826 

Cooray A., 2005, MNRAS, 364, 303 

Cooray A., 2006, MNRAS, 365, 842 

Croton, D.J., et al., 2006, MNRAS,365 ,11 

De Lucia C, et al., 2004, MNRAS, 348, 333 

Diemand J., Kuhlen M., Madau P., 2007, ApJ, 667, 859 

Ferguson A.M.N. , Irwin M.J., Ibata R.A., Lewis G.F., Tanvir N.R., 

2002, AJ, 124, 1452 
Font A.S., Johnston K.V., Bullock J.S., Robertson B.E., 2006, ApJ, 

638, 585 

Gallagher III J.S., Ostriker J. P., 1972, AJ, 77, 288 

Gao L., White S.D.M., Jenkins A., Stoehr F., Springel V., 2004, 

MNRAS, 355, 819 
Giocoli C, Tormen G., van den Bosch F.C., 2008, MNRAS, 386, 

2135 

Gonzalez A.H., Zaritsky D., ZabludofT A.I., 2005, ApJ, 618, 195 
Gonzalez A.H., Zaritsky D., ZabludofT A.I., 2007, ApJ, 666, 147 
Helmi A., White S.D.M., de Zeeuw P.T., Zhao H., 1999, Nature, 
402, 53 

Henriques B.M., Bertone S., Thomas P.A., 2008, MNRAS, 383, 
1649 

Ibata R.A., Gilmore G., Irwin M.J., 199, Nature, 370, 194 
Jing Y.P., Mo H.J., Borner G., 1998, ApJ, 494, 1 
Johnston K.V., Hernquist L., Bolte M., 1996, ApJ, 465, 278 
Johnston K.V., Sackett P.D., Bullock J.S., 2001, ApJ, 557, 137 
Kang X., van den Bosch F.C., 2008, ApJ, 676, LlOl 
Kim J., Park C, Choi Y.Y., 2008, ApJ, 683, 123 
Lin Y.T., Mohr J. J., 2004, ApJ, 617, 879 

Madau P., Ferguson H.C., Dickinson M.E., Giavalisco M., Steidel 

C.C., Fruchter A., 1996, MNRAS, 283, 1388 
Mandelbaum R., Seljak U., KaufTmann G., Hirata CM., 

Brinkmann J., 2006, MNRAS, 368, 715 
Merritt D., 1983, ApJ, 264, 24 

Mihos J.C., 2004, in Clusters of Galaxies: Probes of Cosmological 
Structure and Galaxy Evolution, ed. J. S. Mulchaey, A. Dressier, 
& A. Oemler, 277 



is a significant fraction of ~ 19% (upper-limit) or ^ 5% 
(lower- limit). The implication of the existence of such 
a stellar component for galaxy formation has yet to be 
explored. Indeed, if the total amount of stars in groups, 
clusters and in the universe is larger than that implied 
by the observed stellar mass function of galaxies, and if 
semi-analytical models and numerical simulations do not 
resolve the ICS component, then it would be incorrect 
to use the observed stellar mass/luminosity function of 
galaxies to constrain star formation efficiency and feed- 
back. Furthermore, the fraction of the ICS component 
in halos of different masses may also convey important 
information about the evolution of galaxies in different 
environments. Clearly, more theoretical work is required 
in order to make full use of the information provided by 
the ICS component. 
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